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Abstract. Bacteriophages, phages for short, are viruses of bacteria. The majority 
of phages contain a double-stranded DNA genome packaged in a capsid at a density of 
^ 500 mg/ml. This high density requires substantial compression of the normal B form 
helix, leading to the conjecture that DNA in mature phage virions is under significant 
pressure, and that pressure is used to eject the DNA during infection. A large number 
of theoretical, computer simulation and in vitro experimental studies surrounding this 
conjecture has revealed many — though often isolated and/or contradictory — aspects 
of packaged DNA. This prompts us to present a unified view of the statistical physics 
and thermodynamics of DNA packaged in phage capsids. We argue that the DNA 
in a mature phage is in a (meta)stable state, wherein electrostatic self-repulsion is 
balanced by curvature stress due to confinement in the capsid. We show that in 
addition to the osmotic pressure associated with the packaged DNA and its counterions, 
there are four different pressures within the capsid: pressure on the DNA, hydrostatic 
pressure, the pressure experienced by the capsid, and the pressure associated with the 
chemical potential of DNA ejection. Significantly, we analyze the mechanism of force 
transmission in the packaged DNA, and demonstrate that the pressure on DNA is 
not important for ejection. We derive equations showing a strong hydrostatic pressure 
difference across the capsid shell. We propose that when a phage is triggered to eject by 
interaction with its receptor in vitro, the (thermodynamic) incentive of water molecules 
to enter the phage capsid flushes the DNA out of the capsid. In vivo, the difference 
between the osmotic pressures in the bacterial cell cytoplasm and the culture medium 
similarly results in a water flow that drags the DNA out of the capsid and into the 
bacterial cell. 
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1. Introduction 

Bacteriophages, or pliages, are viruses of bacteria. Phages consist of a protein capsid 
that encapsidates their genome, and a tail — a hollow tube connected to the capsid 
via a portal complex. During infection, the phage tail attaches to a host bacterium, 
punctures the cytoplasmic membrane and its genome translocates through the portal 
and the tail into the bacterial cytoplasm. Infection initiates the phage life-cycle: within 
the bacterial cytoplasm the genome is transcribed and replicated, phage proteins are 
synthesized, and new genome copies are packaged into newly assembled capsids. The 
cycle ends with lysis of the host cell and the release of multiple progeny. Understanding 
the mechanism(s) of phage genome ejection is important, not only in the insights it 
provides on DNA structure but also to provide a model for how eukaryotic viruses may 
release their nucleic acid in to the cytoplasm or nucleus of an infected cell. 

The genome of a mature phage virion is usually a B-form, double-stranded DNA 
(dsDNA). A common phage may have 35-50 kb DNA (DNA and genome will be used 
interchangeably) packed into a capsid of ~ 60 nm diameter, i.e., the phage DNA is 
packaged tightly within the capsid in a condensed state, at a linear compression factor 
~ 250, or at a density ~ 500 mg/ml [T]-[T ^[20l[2T| . More than five decades ago, this 
compression led to the conjecture (henceforth referred to as the "pressure-conjecture") 
that DNA in a mature phage capsid is under significant pressure, and is kept in place 
by means of a "plug": protein(s) in the tail tube. When the plug is opened by the 
action of the appropriate receptor, the pressure of the DNA causes its release into the 
cytoplasm of an infected cell in a biologically passive manner [22| - [25] . The difficulty of 
obtaining experimental kinetic data for the ejection of any phage genome into an infected 
cell resulted in the conjecture becoming "fact". The few published examples that were 
inconsistent with the theory elicited little response and most textbooks simply refer to 
the "DNA injection" step without further elaboration. However, at the beginning of this 
millennium, the pressure conjecture resurfaced in the biophysics community. A single 
molecule study reported that the packaging motor of phage 029 is capable of packaging 
DNA into a phage capsid against a force of magnitude > 60 pN [26]. The force, which 
was assumed to be totally conserved, and thus available for the subsequent ejection step, 
and which therefore can be termed as the "ejection force", was simply defined as the 
pressure on the DNA within the capsid, multiplied by the cross-sectional area of the unit 
cell corresponding to the hexagonal lattice arrangement of the packaged DNA. It was 
estimated that the pressure on the packaged DNA, as well as the pressure on the inside 
of the capsid is in the order of 60 atm for a mature phage virion. This work was quickly 
followed by a large number of computer simulations [271 - I30] . theoretical analyses [3TH35] . 
and in vitro experimental studies [36] - I^H3] - H8] surrounding the pressure-conjecture, a 
significant fraction of these studies stating the 60 atm pressure on the packaged DNA 
for a mature virion as a matter of fact. 

Many recent theoretical, computer simulation and in vitro experimental studies 
have been principally directed to quantify how the pressure on the DNA can be 
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understood from a thermodynamic perspective. However, to the best of our knowledge, 
the issue of the mechanical transmission of the pressure along the packaged DNA helix 
has received no attention. Putting aside that question for the moment, a general 
consensus has been reached on a thermodynamic description of the pressure on the 
DNA. DNA is a charged polymer with persistence length ~ 50 nm; DNA confined 
inside a capsid of diameter ~ 50 nm has a large free energy cost F, relative to the 
state of the DNA outside the capsid. The pressure on the packaged DNA can be 
thermodynamically derived from this free energy. However, fundamental disagreements 
remain between different theories on the thermodynamic origin of the pressure acting on 
the packaged DNA. In light of these disagreements, together with the question of how 
mechanical transmission of force along a flexible polymer, naturally raises the question: 
"Does the pressure-conjecture necessitate a biologically passive ejection force"? 

In order to appreciate the outstanding issues more easily, we provide a brief overview 
of the structure of packaged phage DNA, and existing approaches to describe phage DNA 
packaging and ejection. 

1.1. Structure of the packaged DNA in mature phage virions 

Experiments on the structure of the packaged DNA within phage capsids have a long 
history: preliminary evidence of hexagonal packaging of the DNA for T2 and T7 phages 
was first obtained in 1961 [1], and X-ray diffraction was also used to establish that DNA 
forms concentric layers within the capsid [2H1]. Several theoretical models have been 
proposed for the topology of the packaged DNA: wound into a spool |19H5T], a liquid 
crystal with hairpins |52ll53] or defects [5l], to name but a few, subjecting the topic 
to much debate. The first unambiguous high-quality images of DNA organization in 
mature phage capsids were obtained using cryo-electron microscopy in 1997. Averaged 
over many mature T7 tail-deletion mutant virions, viewed along the portal axis, the 
packaged DNA showed patterns of circular striations, spaced ~ 2.5 nm apart (Fig. 
H]) A computer-modeled projection (side views) of spooled DNA within the capsid, 
again averaged over many virions, strongly suggested that T7 DNA is wrapped axially 
in concentric shells of a toroid. In each shell, DNA is coiled with an axial rise of ~ 2.5 
nm per turn, and the diameter of each turn is imposed by the previous shell. Following 
this pioneering study, a succession of cryo-electron microscopy reconstructions have 
provided support for a toroidal structure: A [5], both isometric and prolate T4 |6l[7], 
P22, el5, 029 [8HI3], KIE and Kl-5 [H], P-SPP7 [15] and N4 [20] all reveal consecutive 
layers of toroidal DNA spaced by ~ 2.5 nm. Together, these observations suggest that, 
when averaged over many virions, the well-defined hexagonal lattice for the packaged 
DNA may be a generic feature of all phages (Fig. [2]). Separately, it should be noted 
that, where it was resolved in these reconstructions, the leading end (first end to enter 
the infected cell) is seen to extend into the narrow portal channel. Furthermore, in 
earlier studies, the leading end of the genome could be cross-linked to the tail of several 
mature phage virions [T614T9] . Portal/tail insertion of the leading DNA end during virion 
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morphogenesis is also likely a general feature as it ensures that the subsequent DNA 
ejection step is efficient. 




Figure 1. Source: Ref. Hj. (b) Cryo-electron images of mature T7 tail-deletion mutant 
virions, viewed along the portal axis; bar corresponds to 25 nm. (c) Image obtained 
by averaging over 21 virions, the closed triangle marks the discontinuity between the 
second and the third DNA-associated rings, (d) Averaging over 77 virions; on the right 
appears the azimuthally averaged image: S is a dense ring corresponding to the capsid, 
and 9 other peaks corresponding to 9 DNA rings. Reproduced with permission from 
Elsevier Inc. 

It is of paramount importance to emphasize that the hexagonal toroidal spool 
structure of the packaged DNA is an average property. This point is illustrated 
beautifully in a recent paper on phage T5 [55]. Further, to be packed within a capsid 
the DNA helix has to cross itself and therefore cannot lie on a perfect lattice. The 
toroidal structure has also been suggested to hold only for the part of the DNA close to 
the inner wall of the capsid, and not for the entire DNA (e.g., Refs. |21]|29]). 




Figure 2. (Left) Source: Ref. [32]. Idealized DNA packaging in a coaxial toroidal form 
within the capsid. Reproduced with permission from Elsevier Inc. (Right) Hexagonal 
lattice formation by the DNA spool, within the capsid. The black filled circles at the 
center of each hexagonal unit cell (an example is painted pinkish red) denotes the 
cross-section of the DNA. 

1.2. Existing approaches for the thermodynamics of packaging and ejection forces 

There have been two main — and mutually conflicting — approaches that address 
the thermodynamic origin of packaging and ejection forces: continuum mechanics and 
coarse-grained molecular mechanics models. In the continuum mechanics model [3TH35] , 
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packaged DNA is seen as a uniformly charged rod, with a certain persistent length Ip, 
organized in a toroidal spool on a hexagonal lattice (Fig. |2]). The continuum mechanics 
model is simple enough to allow analytical calculations: packaged DNA is assumed to 
repulsively self- interact (with electrostatic Coulomb interaction), allowing calculation of 
the stored electrostatic energy Up^^\ Similarly, the bending energy U^'' is obtained from 
the toroidal spool arrangement of DNA on a hexagonal lattice. Both the electrostatic 
and bending energies are assumed to be zero for free DNA outside the capsid. Thereafter, 
with the assumption that the difference between the entropy of free and packaged DNA 
is negligible, the difference between the Helmholtz free energy of packaged and free DNA 
is obtained as F = uj^^^ + Note that in a given buffered solution and a capsid of 
volume V, made of a perfectly rigid material, F is only a function of the length L of the 
DNA within the capsid, i.e., F = F{L). In particular, F increases as a function of L, 
and it is —dF/{Aj)^AdL), the pressure associated with the chemical potential of DNA 
ejection (the chemical potential of DNA ejection is simply defined as the negative rate 
of change of the free energy F with respect to the length of DNA within the capsid) 
that manifests as the ejection pressure, where A^na is the cross sectional area of the 
DNA. Note that OF/dL has the units of force; consequently, in the continuum mechanics 
model, —dF/dL is viewed as the ejection force, and OF/dL as the packaging force (note 
that —OF/dL is also the chemical potential of DNA ejection, as defined above). 

In the in vitro A experiments [361439] . phage virions are immersed in a solution 
containing PEG and/or DNA condensing agents, and DNA is ejected when triggered by 
the LamB receptor protein. The length of the DNA that remains within the capsid at 
the end of the ejection process depends on the concentration of the PEG and/or DNA 
condensing agent. (The basic principles underlying these in vitro experiments were 
first presented in a theoretical paper more than four decades ago [25]). The continuum 
mechanics model has been widely advertised to explain the in vitro ejection data, not 
just for A but for all phages — both in vitro and in vivo. A somewhat deeper look, 
however, reveals that the agreement between theory and the in vitro experiments using 
T5 (one of only three phages studied in osmotic suppression experiments), remains 
poor [iniHT]. The model also neglects how the ejection of proteins from infecting virions 
(which is a feature of most phages) into the cell is accomplished, and it cannot therefore 
describe the complete infection process. Further, there is also disagreement between the 
continuum mechanics theory and the experimental kinetics of T7 DNA translocation 
into the bacterial cytoplasm in vivo [56l|57]. All but the leading 0.5 kb of the 70 kb 
genome of the unrelated phage N4 is also known to be internalized by transcription 
in the cell [20l|42], and marker rescue experiments with phage SP82 suggest that the 
phage genome enters a cell at a constant rate [13]. None of the in vivo experiments 
that show kinetics of genome internalization are consistent with theories derived from 
in vitro ejection studies. 

The continuum mechanics model suffers from fundamentally unrealistic assump- 
tions that diminish its value. In particular, free energy is an equilibrium concept] ob- 
taining an ejection pressure by —dF/{A-Q^p^dL), in order to explain the experimental 
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data on the irreversible ejection process, has the underlying assumption that DNA re- 
mains at equilibrium in a toroidal spool configuration on a hexagonal lattice at all stages 
of ejection. Not only is this assumption unsubstantiated by polymer physics theory, it 
is also inconsistent with actual experimental data on T5 [55j ; the latter is illustrated in 
Fig. [3] below. In order to circumvent some of these problems, the continuum mechanics 
model invokes two free parameters to fit experimental data. 




Figure 3. Cryo-EM images of a single T5 virion at different stages of ejection: (a) 
full capsid to (h) empty capsid, showing the lack of toroidal spool configuration of the 
DNA hexagonal lattice at all stages of ejection, in contrast to the assumption made by 
the continuum mechanics model. Figure reproduced from Ref. |55) , with permission 
from Elsevier Inc. 

In coarse-grained molecular mechanics models [271430] . DNA is modeled as a 
polymer, with a persistence length Ip, and polymer dynamics is allowed to take its 
own course in simulations of phage genome packaging and ejection. The model does 
not allow analytical calculations — it is effected only by computer simulations, and 
in work published to date, the DNA is modeled without surrounding water and 
counterions. However, there is no assumption that DNA remains at equilibrium in 
a toroidal spool configuration on a hexagonal lattice during packaging or ejection. The 
main limitation of coarse-grained molecular mechanics model is that because processes 
are only materialized by computer simulation, they are difficult to correlate with 
experimental data. 

1.3. This paper 

In summary, our understanding of the thermodynamics of phage DNA packaging and 
ejection, provided by the existing literature, is far from satisfactory. In this paper 
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we present a comprehensive treatise of the thermodynamics of packaged DNA and the 
mechanism of genome ejection, independent of any specific model. In Sec. [2] we show 
that there are four distinct pressures within the capsid: pressure on DNA, hydrostatic 
pressure, pressure experienced by the capsid, and pressure associated with the chemical 
potential of DNA ejection. These are all different thermodynamic quantities. We 
also show that for amature phage the hydrostatic pressure within the capsid is much 
higher than outside, and that the pressure on the DNA is a lot smaller than commonly 
envisaged. In Sec. [3] we take up the issue of force transmission along the DNA helix, 
and show that the pressure on DNA is not important for genome ejection. Instead, in 
Sec. m we suggest that when a phage is triggered in vitro to eject by interaction with its 
receptor, the (thermodynamic) incentive for water to enter the capsid flushes the DNA 
out. We argue that in vivo, the difference between the osmotic pressure within the cell 
cytoplasm and the outside culture medium initiates, and importantly maintains, a flow 
of water from the culture medium into the capsid, and then from the capsid down the 
phage tail-tube into the bacterial cell cytoplasm, dragging the DNA out of the capsid 
and into the cell. We show that this theory is consistent with both studies of in vitro 
ejection and observations of infection in vivo. 



2. Thermodynamics of phages in vitro 
2.1. Entropy of DNA ejection in vitro 

We consider an experiment with an ensemble of realizations, each containing a phage 
capsid initially placed in a certain buffer solution: for each realization a DNA of length 
L is packaged within a capsid, whose volume we denote by V . The DNA is allowed to 
exit the phage in vitro, and at the end of the experiment the entire genome ends up in 
the buffer. The experimental system is kept at a temperature T at all times. The only 
interaction of the system with the environment is to exchange (thermal) energy. For the 
system in its entirety, we define the initial (DNA within the capsid) free-energy and 
the entropy S^; the final (DNA in the buffer) free-energy Funp and entropy Sunp. The 
quantities Fejection = -^unp - -^p, and ^ejection = 5'unp - Sp are then the free energy and 
entropy of ejection, respectively. 

The above isothermal experiment, when actually performed in a laboratory, is 
exothermic [58j. Earlier experiments using differential scanning calorimetry led to the 
same conclusion [59146T] : i.e., the system releases an amount of heat, of magnitude 
Rejection to the environment. This tells us that 5'ejection < 0; however, the amount of heat 
released to the environment cannot give us the value 5'ejection because the ejection process 
is spontaneous and irreversible. All we can say, using the second law of thermodynamics, 

is that 1 5'ejection I ^ Qejection/^- 

The result that the entropy of ejection is negative, from the point of view of DNA 
confinement in the capsid, is counter-intuitive. One would expect the configurational 
entropy of the DNA to be strongly reduced due to confinement when compared to the 
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same DNA in solution; in other words, from the configuration entropy of the DNA 
in isolation one would expect S'ejection > 0. This is, in fact, correct; however, the DNA 
cannot be considered in isolation. The configurational entropy of the DNA is only a very 
small part of the entropy of the entire system; most of the change in entropy associated 
with ejection of DNA from its confined state within the capsid into the environment 
involves the ordering of water molecules around the free DNA [62]. DNA, being a 
charged molecule, orders the dipole orientation of nearby water molecules. There are 
many fewer water molecules surrounding the DNA inside the capsid when compared 
to the same DNA in solution. Consequently, although the DNA molecule does gain 
entropy following ejection from the capsid, so many more water molecules lose entropy 
that the entropy of the entire system decreases as DNA is ejected from the capsid. 

This issue can be used to evaluate the current approaches for analyzing the physics 
of phage ejection. DNA configurational entropy is not considered a significant parameter 
in the continuum mechanics model, while molecular mechanics models [271430] show 
that the configurational entropy constitutes a sizable fraction of the free energy of the 
packaged DNA. 

2.2. Thermodynamics for the capsid content of a mature phage in a buffer solution 

Most phage capsid shells are fully permeable to water and small ions or molecules, but 
are impermeable to large molecules in the buffer or to the encapsidated DNA. Some 
phage capsids, like those of T4 and its relatives, are much less permeable, even to small 
molecules or ions, including cesium chloride and ammonium acetate. Such capsids can 
easily be broken by osmotic shock, a phenomenon that resulted in the classic electron 
micrograph showing phage T2 DNA outside a ruptured virion [63] . 

We assume that the entire system, with DNA packaged within the capsid, can 
be thermodynamically separated into two subsystems separated by the capsid shell. 
The subsystem within the capsid consists of DNA (perhaps in addition to protein 
molecules, which we will not specifically refer to further), and an aqueous environment 
containing small solutes (including ions) that can permeate the capsid shell, while the 
subsystem outside usually consists of water and solute molecules that can or cannot 
permeate the capsid shell. Apart from exchanges of solutes that can permeate the 
capsid shell, these two subsystems do not interact. This assumption allows us to express 
Fp = -Fp^in + Fp^out, where Fp^in and Fp^out are, respectively, the free energy of the contents 
of the capsid and that of the environment outside the capsid. Clearly, Fp in is a function 
of the equilibrium capsid volume V, DNA length L, number of water molecules {riw.in} 
and the number of solute molecules of all species {ris.m} within the capsid. Similarly, 
Fpout is a function of number of water molecules {n^^out} and the number of solute 
molecules of all species {ns,out} in the solution. Both subsystems are "open" in the 
thermodynamic sense, since they exchange water and small solute molecules. We assume 
(realistically) that the volume of the buffer in which the phage is immersed is ^ V, so for 
analyzing the thermodynamics of the phage and its contents, the external solution can 
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be considered to be an infinite reservoir for water and small permeable solute molecules. 
The number of solute molecules within the phage capsid {us^in} is then fixed by the 
chemical potential of the small solute molecules in the buffer, collectively denoted by 
{/is} (for the case of counterions, via the Donnan equilibrium). This simply means 
that Fp^in = Fp^in(V, The free energy Fp^in(V, L)|{^g} can be further dissociated 

into the interaction free energy F^™*^(V, L)|{^g} of the ions, DNA and water molecules 
within the capsid, and the DNA configurational free energy (this includes bending, or 
curvature, energy and the configurational entropy of the DNA) Fc{V, L)\^^^y of the DNA 
is then given by 



We also assume that water is incompressible and that small solutes do not occupy any 
physical volume. With these assumptions, in equilibrium, the number of water molecules 
within the capsid is then determined by V minus the physical volume of the DNA, and 
we can now explicitly focus on a number of thermodynamic issues. 

2.2.1. Pressure associated with the chemical potential for DNA ejection If we assume 
that the capsid is made of a perfectly rigid material, then from the above definitions 
Fp_in(V, L)|{/i4 = Fp^in(-^)|{/i4, we can define the pressure associated with the chemical 
potential for DNA ejection. Consider two situations using the same buffer conditions 
(i.e., at fixed chemical potentials of permeable ions), one where DNA of length L and 
the other a length of L — dL, is packaged within the capsid. We assume that Donnan 
equilibrium conditions are satisfied for both lengths of DNA within the phage capsid. 
One can then define the pressure associated with the chemical potential of DNA ejection 
by comparing the free energy of the capsid content under these two situations, viz. 
^P,m(^)|{M4 and Fp,in(^ - c^^) I as: 



where Adna is the cross-sectional area of the DNA. The pressure Pejection is the pressure 
that the plug at the end of the tail tube feels from within the phage. 

2.3. Pressure on the capsid from within, hydrostatic pressure imbalance across the 
capsid and pressure on the DNA 

In Sec. U.2\ we considered the thermodynamics of the capsid content as a subsystem, 
under a fixed chemical potential of permeable small solute molecules in the buffer, while 
assuming that the number of water molecules within the capsid is simply determined by 
the capsid volume. This is however an incomplete description of the thermodynamics 
of the entire system (namely the phage immersed in the buffer solution), as we did 
not consider the thermodjTiamic penalties associated with exchange of water molecules 
across the capsid shell. If the capsid is made of a perfectly rigid material without 
the possibility of expansion or contraction, such a description would be correct, but in 



Fp,i,(r,L)|{^3l = Ff *)(r,L)|{^3} + F,(V^,L)||^3|. 



(1) 



ejection 




(2) 
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reality, the capsid is comprised of protein molecules, and therefore has a finite, albeit 
high, rigidity [65]. Keeping the assumption that water is incompressible, rather than 
a fixed capsid volume determining the number of water molecules internalized, the 
capsid volume should actually be determined by the thermodynamics of water molecule 
exchange across the capsid shell, which is equivalent to a semi-permeable membrane. 

We now argue that for a given length L of the DNA within the capsid, the 
equilibrium volume V is determined by the buffer composition. We presuppose that 
an expansion/contraction of the capsid allows for an exchange of a small number of 
water molecules, relative to the total number of water molecules present in the external 
buffer. We can then safely assume that the chemical potential of permeable small solutes 
in the buffer remains unaltered by the exchange of water molecules across the capsid 
shell. 

First, the pressure on the capsid from within, Pcapsid, can be calculated from the 
V -dependence of the free energy Fp in(V, as follows. Imagine a virtual, uniform 

expansion of the capsid leading to a volume increase from V toV + dV (the extra volume 
will be filled by water from the buffer and permeable solute molecules as dictated by 
{/is} — for the case of ions, via the Donnan equilibrium). For a given length of DNA, 
under fixed {/is}, this leads to a new value of the free energy of the inside material of 
the capsid, -Fp,m(^ + dV, but also, since the buffer solution loses a volume dV of 

water, the free energy of the buffer solution increases by an amount Tio^tdV , where vTout 
is the osmotic pressure of the buffer solution. The pressure on the capsid from within, 
is therefore given by 

-'capsid "out' \^) 

Secondly, water's thermodynamic incentive to enter the capsid is given by the 
osmotic gradient across the capsid shell, namely [vTin — vTout], where vTin is the osmotic 
pressure within the capsid shell, and is defined by [see Eq. ([1])] 

dv ■ ^> 

Simultaneously, when entering the capsid, the water molecules exit a zone of hydrostatic 
pressure Phydm.out, the hydrostatic pressure in the solution outside the capsid, and enter 
the capsid which may have, in principle, a different pressure, Phydro,in- The equilibrium 
capsid volume is then determined by the condition that [vTin — vTout] is counter-balanced 
by the pressure-volume work done by water molecules in entering the capsid 

-Phydro,in -Phydro,out ^in ^out- (5) 

It is interesting to note that Eq. (jS]) is well-known to hold for bacterial and plant cells: 
the osmotic pressure gradient is counter-balanced by a hydrostatic pressure differential 
(turgor) [68]. Turgor allows cells to enlarge and thus facilitates growth. 

Further, an appreciation for Eq. ([5]) is afforded by a straightforward gedanken 
experiment that connects to elementary physical chemistry. Consider a vertical U- 
tube with a membrane at the lowest point separating the two arms. The membrane is 
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permeable to water, but impermeable to, say, sugar molecules. We fill up the two arms 
of the U-tube to equal height: the left arm with pure water, and the right arm with 
sugar solution. As time progresses, water from the left arm will permeate into the right 
arm, reducing the height of the water column on the left, while increasing the height 
of the sugar solution column on the right. At equilibrium, the height of water in the 
left column will be lower than the sugar solution, meaning that there is a hydrostatic 
pressure difference P across the membrane. Moreover, at equilibrium the concentration 
of sugar in the right column is non-zero, while it remains zero in the left column by 
construction. There is therefore an osmotic gradient across the membrane, which is 
precisely counter-balanced by the hydrostatic pressure P. Equation (jS]) describes the 
same equilibrium, where the difference between the osmotic pressures inside and outside 
the capsid is counter-balanced by the hydrostatic pressure difference across its protein 
shell. 

Equations ([3]) and ([5]) allow us to appreciate not only how the pressure on the capsid 
of finite rigidity is mechanically materialized, but also how the pressure Pdna that the 
capsid mechanically transmits to the DNA (and by action-reaction, the pressure the 
DNA transmits to the capsid). First, the force balance on a small capsid surface tells 
us that the capsid pressure Pcapsid equals the hydrostatic pressure gradient across the 
capsid plus the pressure that is mechanically transmitted to the DNA inside; i.e.. 



Note in Eq. ([7]) that the osmotic pressure within the capsid tT[^ is derived from the 
interaction energy of the ions, DNA and water molecules within the capsid, while the 
same interaction energy appears in Fp^in(V, but with opposite sign [Eq. (jlj)]. 

Consequently, the net contribution of this interaction energy to the r.h.s. of Eq. ([7]) is 
zero. This conclusion implies that the pressure on the DNA mechanically transmitted 
by the capsid comes entirely from the stiffness of the DNA. We verify this conclusion 
in Sec. [3] for a single turn of the DNA. Therein we also discuss why Eq. ([7]) is not in 
conflict with the experiment by Smith et al. [26]. 

It is straightforward to show that Eqs. (jSjlZl) are consistent with the Gibbs free 
energy minimization of the entire system (the mature phage plus the buffer solution). 
Moreover, we note that due to the cancellation of the interaction energy term in Eq. (7) 
as explained above, the compressive forces on the DNA mechanically transmitted by the 
capsid are small compared to Pcapsid^ the bulk of Pcapsid is mechanically transmitted to 
the water inside. This is in fact confirmed from the structural and interaction properties 
of both the packaged DNA and the structure of the protein molecules comprising the 
capsid shell. Raman spectral studies of mature P22 and T7 phages have not found any 
evidence of any structural alteration of protein or DNA that would be expected if they 
were subjected to high pressure [66l|67]. Only the configuration of the phosphodiester 



-f capsid -f hydro, in -Phydro,out ~l~ -PdNA' 

Incorporating Eqs. ([1]) and ([1116]), this reads 



(6) 



_ gFp,,(V,L)||^4,^ 
Pdna - 



dV 



(7) 
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groups in packaged DNA are perturbed from that found in free DNA solutions, and the 
capsid protein structures of mature virions are indistinguishable from empty, DNA-free, 
particles. 

Finally, we note from Eqs. ([2117]) that the osmotic pressure vTin within the capsid, 
the pressure Pejection associated with the chemical potential of DNA ejection, and the 
pressure Pdna imparted by the capsid on the DNA are entirely different thermodynamic 
quantities. Importantly, it must be appreciated that these pressures are all equilibrium 
quantities, and they cannot be used to describe a non-equilibrium situation. 

3. Can 

^DNA and Pejection plciy any role in DNA ejection dynamics? 

In this section, we investigate whether Pdna and Pejection play a role in DNA ejection. 
We first discuss the case of Pdna within the context of the continuum mechanics model; 
in order to do so we start by analyzing the stability of the DNA within a phage capsid. 

3.1. Stability of the DNA within a phage capsid 

As described in Sees. 11.21 continuum mechanics model starts with the assumption that 
the DNA, within a capsid made of a perfectly rigid material, takes a toroidal spool 
configuration on a hexagonal lattice at all stages of ejection. Given that the outer 
radius of the spool is fixed by the capsid dimensions, all structural aspects of the spool 
are then determined by the spacing of the hexagonal lattice. Phages with internal cores 
constitute a special case: the radius of the innermost ring of the spool is simply given 
by the radius of the core. For a given ionic condition {rig in} inside the capsid — which 
in turn is fixed by the chemical potential {/Xg} of the small solutes in the external buffer 
(for the case of counterions, via the Donnan equilibrium) — and DNA of length L within 
a capsid of volume V, the free energy of the material inside the capsid is a function of 
an alone, where an is the hexagonal lattice spacing; i.e., Pp,in(V, = F{aH)\{fis}- 

Note that the configurational entropy of the DNA is not a meaningful quantity in the 
continuum mechanics model, hence in this section Fc only contains curvature (bending) 
energy of the DNA; i.e., F^{V, L)\{f,^y = U.^an). 

Within the continuum mechanics model, one can argue that an is fixed, via a trade- 
off between P'^™*)(a//)|{^g} and Ucian) [see Eq. (IT])]. The interaction free energy can 
be lowered by increasing an- However, since the "volume constraint" (meaning that 
the total available volume within the capsid is V) must be obeyed, increasing an will 
reduce the minimum radius of the DNA spool, because within the inner part of the 
spool the DNA will be much more tightly bent, which will increase its bending energy. 
In this manner, a reduction in interaction free energy increases curvature energy and 
vice versa. In other words the trade-off between these two quantities dictates that: 

^ ^ — ^ ~ y^) 

which minimizes the free energy P(a//) ||^^}. We denote the value of qh-, obtained from 
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Eq. (IH]) by a^'' , and the corresponding minimum of the free energy of the inside content 



of the capsid by -Fp^inliAtj = {fMs}- If it is further assumed that most of the 

interaction free energy variation due to variations in caused by a virtual uniform 
expansion of the capsid [as in Eq. (^], involves the DNA helix, then the osmotic pressure 
within the capsid can also be defined as 



[a 



TTi, 



H JI{m4 



3^/3a\j L da 



(0) 
H 



(9) 



where the volume that the DNA toroid occupies within the hexagonal lattice is given 



by 3x/3 



'■H 



L/4. This analysis for the specific case of the phage T7 — a phage with 



an internal proteinaceous core — can be found in Ref. [& 



3.2. Can Pdna transmit down the tail tube? 

The only way Pqna can play a role in DNA ejection is if it can transmit from within 
the capsid, where the bulk of the DNA is located, down the tail tube. To evaluate this 
idea, we begin with the appreciation that Pdna favors a uniform expansion of the capsid. 
However, since DNA is not a continuum material, Pdna does not trivially transmit along 
the helical axis into the tail tube. Nevertheless, when the DNA organization within the 
capsid is averaged over many phage particles, the picture that emerges is that pressure 
on the toroidal spool of DNA does translate into a compressive force / acting along the 
contour of the DNA. Following Sec. 13. H we start by analyzing this force within the 
continuum mechanics model, and we then ask whether this force can transmit down the 
tail tube. 



(a) 




pitch of the DNA toroid 




Figure 4. (a) One full turn of the DNA toroid around the axis of the idealized toroidal 
spool in Fig. [2] (b) With the pitch ignored: the corresponding DNA ring of radius 
R, and the compressive force / along the contour of the DNA, shown on a differential 
length element Al of the ring. 



We consider one full turn of DNA around the axis of the idealized toroidal spool 
[Fig. Ht^a)]. For simplicity, we ignore the pitch of this turn, i.e., it is idealized as a DNA 
ring of radius R, and we focus on a differential length segment Al. The DNA ring and 
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the compressive force along the DNA contour are shown in Fig. ID^b). The fact that 
the entire DNA spool is under a pressure Pdna within the phage capsid implies that 
the isolated ring itself is under a pressure P, but whose relation to Pdna we have not 
specified. 

Within the scheme of continuum mechanics, the force balance equation for this 
differential length segment along the transverse direction (along a line joining the 
differential length element and the center of the ring) is given by 

f sine ^ fe = PAH, (10) 

where t is the cross- sectional diameter of the DNA. Further, the geometrical con- 
sideration from Fig. Hl^b) tells us that 

sine ^ e = Al/R, (11) 

i.e., 

f = PRt. (12) 

We now note that, just like in Eq. ([8]), the DNA ring stays in its configuration by 
the balance of two competing thermodynamic forces, one derived from its curvature 
energy Uc, the other from the energy of self-interaction with nearby DNA strands. The 
thermodynamic force related to curvature energy favors an increase of the radius of the 
ring, while the thermodynamic force related to the self-interaction free energy favors a 
decrease. The magnitude of this force is given by 

(OrrmPi 27l{kBT)lp 

implying that 

f = PRt=^-^. (14) 

At room temperature (T = 300 K) ksT = 4.14 x 10~^^ Joules. If we take Ip = 50 nm, 
the widely used value, then 

, 2.07 X 10-28 2 

/ = Nm2. (15) 

Thus, depending on the radius of the DNA ring, the magnitude of the compressive force 
/ along the contour of the DNA can be a few pN. Equation (fTSll tells us that the pressure 
on this DNA ring is only in terms of the stiffness of the DNA, as we argued in Sec. 12.31 
How much of this force can transmit down the tail tube? We should remember that 
force is a vector quantity. Although the above exercise shows that there is a compressive 
force along the contour of the DNA ring that make up the toroidal spool, it is necessary 
to follow the actual turns of the DNA leading into the tail tube in order to identify how 
much of the force / is actually being transmitted. However, the definition of persistence 
length Ip in polymer physics means that while a DNA molecule of length < Ip is semi- 
flexible, a DNA molecule with length ^ Ip is simply a flexible polymer. It is therefore 
impossible to transmit compressive forces along the contour of a DNA that is ^ 50 nm. 
This consideration alone ensures that any compressive force along the contour of, e.g.. 
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an ~ 16 /iin long A genome will not transmit into the tail tube. Therefore, compressive 
forces on the packaged phage genomic DNA cannot be important for DNA ejection. 

Our estimate of / having a magnitude of a few pN is, at first blush, in apparent 
contradiction with the experiment by Smith et al. [26]. A deeper look however 
immediately reveals that Eq. (JT5l) and the forces exerted by the packaging motor are 
actually compatible. Thermodynamically, we have identified four different pressures 
within the capsid: pressure on the DNA Pdna, hydrostatic pressure Phydro.in, the pressure 
experienced by the capsid Pcapsid, and the pressure Pejection associated with the chemical 
potential of DNA ejection. These are all different quantities. Smith et al. [26] measure 
-Pejection; they assumed that the work done against Pejection by the packaging motor is 
used to pressurize the DNA and is conserved in the capsid. However, most of the work 
done against Pejection by the packaging motor is expended in increasing the osmotic 
pressure within the capsid. This increase is due to the expulsion of water molecules 
out of the capsid — a reverse osmosis process — in order to allow the phage genome 
being packaged to condense, while only a small part of the work is used to increase 
Pdna The energy expended during water expulsion is not conserved by the packaged 
genome, which explains our lower estimate of Pdna inside the capsid. In this context 
it is important to remember that osmotic pressure is only a measure for the chemical 
potential of water: the fact that the osmotic pressure inside the capsid is high does not 
imply that the DNA is under enhanced pressure. 

3.3. Is Pejection important for DNA ejection dynamics? 

Equation ([2]) shows that for a mature phage in a given buffer, the tail plug feels 
a thermodynamic pressure Pejection from inside the capsid. Given that Pejection is a 
thermodynamically derived quantity (i.e., derivative of Pp^in with respect to L), it is 
clear that if genome ejection were a quasi-equilibrium process, it would indeed govern 
the ejection dynamics. However, the underlying assumption of quasi-equilibrium is that 
at all stages of ejection the DNA remaining within the phage head stays in a toroidal 
loop, with an ever-increasing spacing between the toroidal strands as ejection proceeds. 
This is the major thesis of the continuum mechanics model, which uses Pejection to explain 
the in vitro experimental data using phage A. 

Whether this thesis and its underlying assumption are correct or not lies in the 
question of time scales. If the time-scales associated with the equilibration of DNA 
(in the form of a toroidal spool on a hexagonal lattice) and the small solute molecules 
within the phage capsid (maintained at a fixed chemical potential by the external buffer), 
are much smaller than the ejection time, equilibrium pressures and forces can be used 
to describe genome ejection. However, recent experimental evidence with T5 clearly 
demonstrates that as the DNA leaves the capsid, it goes through a series of phase 
transitions [55]. The cryo-electron microscopic images of T5 virions during DNA ejection 
strongly suggest that the process is non-equilibrium, which invalidates the use of Pejection 
to describe its dynamics. Comparable experiments have not yet been conducted with 
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other phages, but there is no justifiable reason to assume that the genome of other 
phages, in particular those lacking a defined internal core structure, whose DNA is 
packaged to the same density, would behave otherwise. We are not aware of experimental 
studies that probe the equilibration time-scales of confined DNA, such as within a phage 
capsid. However, the high velocity of phage DNA ejection measured in vitro: up to 75 
kb/sec (T5) with long pauses between distinct steps [Hj, and up to 60 kb/sec (A) in 10 
mM Na"*" buffer or ~ 20 kb/sec in 10 mM Mg^"*" [15], makes it highly unlikely that these 
are quasi-equilibrium processes. These experimental data suggest that Pejection cannot 
be used to describe DNA ejection dynamics even in vitro. 

4. An alternative mechanism of DNA ejection from phage virions 

Sections [2] and [3] illustrate the fundamental problems with the idea that the thermo- 
dynamically derived ejection pressure Pejection (or the corresponding ejection force 
^DNA-Pejection) causcs DNA cjcctiou from phages. There is therefore a clear need for 
a mechanism that can explain the physics of DNA ejection, one that is consistent with 
both in vitro and in vivo experimental data. We suggest such a mechanism below. 

4.I. DNA ejection in vitro 

Of the two main approaches for DNA ejection: continuum mechanics and coarse-grained 
molecular mechanics models, only the former addresses, using Pejection, the DNA ejection 
mechanism in a way that can be compared to in vitro experimental data. In these 
creative experiments, phages are immersed in a solution containing PEG and/or DNA 
condensing agents, and DNA is ejected from the virion when triggered by the receptor 
protein. The length of DNA that remains within the capsid at the end of the ejection 
process depends on the concentration of the PEG and/or DNA condensing agent. 

With the caveats we have discussed concerning the validity of using Pejection 
to describe phage DNA ejection, and remembering that Pejection also involves fitted 
parameters, the continuum mechanics theory does provide good agreement with in 
vitro genome ejection data for A and SPPl (e.g., Refs. [36] - [38| H7 t H8]). However, the 
data obtained with just two phage systems, along with the physics of the continuum 
mechanics model, have often been extrapolated to all phages and in addition to explain 
all phage DNA ejection in vivo. These generalizations ignore data from the third phage 
system that has been studied in vitro, which provides rather different conclusions. At 
the presumed pressures internal to the largely full T5 capsid, in vitro experimental data 
can be fitted to the continuum mechanics model, but at low to moderate pressures, when 
approximately half or less of the genome remain in the capsid, multiple populations are 
found to co-exist |10]: some phages have completely ejected their genomes, whereas 
others have ejected a varying amount of DNA that is not dependent on the external 
osmotic pressure. Furthermore, at < 2 atm external osmotic pressure, most T5 virions 
completely eject their DNA. In contrast, the same external pressure prevents ejection 
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of ~ 40% of the A genome |37]. Theories dependent on Pejection are clearly unable to 
explain how all phages eject their DNA, even in vitro. 



4- 1.1. Hydrodynamic model of phage DNA ejection As discussed in Sec. I3.3t in an in 
vitro experiment when a mature phage is equilibrated in a certain buffer solution, the 
difference between vTin and the osmotic pressure vTout of the buffer solution determines the 
thermodynamic incentive for water molecules to enter the capsid. Equation ([5]) shows 
that if the phage tail is plugged, this incentive is counter-balanced by a hydrostatic 
pressure gradient across the capsid shell (and the capsid, being of high but not infinite 
rigidity, will be slightly expanded from its relaxed state under the enhanced hydrostatic 
pressure from inside). We propose that it is this thermodynamic incentive of water 
molecules to enter the capsid that flushes the DNA out when the tail plug is opened 
by the action of the appropriate receptor. We refer to this idea as the hydrodynamic 
model of DNA ejection. 

When the tail plug is opened by the appropriate receptor, e.g., LamB protein in 
the case of A, FhuA for T5, or YueB780 for SPPl, it is likely that the excess hydrostatic 
pressure within the capsid is partially relieved by a transient water flow down the tail 
tube and, and as a result, the capsid will reduce in size. Although water flow down the 
tail-tube will exert a hydrodynamic force that is likely to drag along the DNA, one end 
of which is already inserted into the tail-tube, this process will not last long enough 
to drag the entire genome out of the virion. While a very small contraction of the 
capsid is unlikely to be detected even by state-of-the-art experiments, we can provide a 
"guesstimate" for how much DNA could be dragged out by this transient water release. 
With a typical capsid diameter ~ 60 nm (assumed spherical), let the capsid diameter 
contract by 0.1 nm; the water released is 47r(30)^(0.05) nm'^. With a typical tail tube 
inner diameter ~ 4 nm, its cross-sectional area is vr(2^) nm^ (ignoring the cross-section 
of the DNA within the tail tube) and the length of the water column passing into the 
tail-tube = 47r(30^)(0.05)/[7r(2^)] nm = 45 nm. Assuming the DNA flows freely with 
this water column down the tail, in this example the DNA, which occupies 25% of the 
cross-sectional area of the tube, can be moved 60 nm, corresponding to ~ 180 base pair 
of a phage genome. 

Once the water pressure difference across the capsid shell is relieved by this transient 
water flow, so long as the difference between the osmotic pressure of the remaining DNA 
within the capsid and vTout stays positive, water will have a thermodynamic incentive to 
move into the capsid using whichever path it can flnd. However, since the capsid has a 
finite volume, any water movement into the capsid can only occur at the expense of an 
equal volume of water, now containing DNA (as a solute) being ejected through the tail 
tube. In other words, following the initial transient water release down the tail-tube, 
continued water movement from the buffer into the capsid up an osmotic gradient will 
cause the DNA to be ejected. This process will stop when the osmotic pressure of the 
leftover amount of DNA in the capsid equals TTout- The physics behind this in vitro DNA 
ejection experiment (omitting the transient water flow down the tail tube due to the 
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excess hydrostatic pressure) is schematically shown in Fig. O 

The hydro dynamic mechanism for DNA ejection commensurate with in vitro 
experimental data. All that is required to drive DNA ejection is an incentive for the 
water molecules to enter the capsid, which is decided dynamically in a non-equilibrium 
manner. There is no a priori assumption that phage genome ejection is a quasi- 
equilibrium process, and thus the model is not dependent on thermodynamic analyses. 



(a) 



phage capsid 



tail tube 




water entering through the capsid shell 
(b) 



in vitro DNA ejection experiment 




flushed out DNA 



tail-plug closed 



tail-plug open 




Figure 5. Schematic of an in vitro DNA ejection experiment, (a) Mature phage 
equilibrated in a certain buffer solution. The difference between TTjn, the pressure on 
the DNA, and the osmotic pressure TTout of the buffer solution is counter-balanced by 
a hydrostatic pressure difference across the capsid shell [Eq. ([5])]. (b) The tail plug is 
opened by the action of the receptor, following any transient water release down the 
tail-tube (not shown), water enters the phage capsid and flushes the DNA into the 
buffer. 

The difference between the osmotic pressure of the residual DNA in the capsid and 
TTout is controlled by two aspects: (i) the presence of any small DNA condensing agents 
that reduces the osmotic pressure of the capsid content without altering vTout, and (ii) 
the presence in the external buffer of large solute molecules like PEG that increase TTout , 
but do not affect the inside osmotic pressure as they cannot penetrate into the capsid. 
Both achieve the same thermodynamic incentive differential for water between inside 
and outside the capsid, and are qualitatively in agreement with in vitro DNA ejection 
data [3&l - [4Tll44] - [48] . When the osmotic pressure of any DNA remaining in the capsid 
= TTout (such as at the end of the ejection process), there will be no further flow of water 
into the capsid, and thus no more DNA ejection. 

The hydrodynamic model requires that the time-scale of water flow into the capsid 
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be sufficiently fast to match the time-scales for DNA ejection in vitro. To this end, 
we refer to an experiment where a mature wild-type A capsid was locally deformed 
(compressed) by up to 25% of its original volume with an AFM tip [65]. Compression of 
the capsid without it rupturing has to drive out internal water and small solutes. When 
the AFM tip was removed, the virion capsid returned to its normal size and shape 
within ~ 4 ms. This suggests that the large surface area to volume ratio of the capsid 
allows copious amounts of water to diffuse — a slow process relative to hydrodynamic 
processes — across the capsid shell in a very short time. Extrapolating these data to 
yield a complete exchange of the contents of the capsid with the external buffer gives 
20 ms; assuming that DNA moves freely, this value results in a rate of phage genome 
ejection in vitro at least 50 times the maximum measured to date. 

By itself, the hydrodynamic model cannot explain either the stepwise ejection of 
DNA from T5 particles or the coexistence of virions containing different amounts of DNA 
at low external pressures. A partial ad hoc solution for T5 (one that is also necessary for 
the Pejection modcl of DNA ejection) is that the tail tube becomes temporarily blocked 
by protein conformational changes, thereby impeding water flow — and thus stopping 
DNA ejection. Note that water flow from the external buffer into the capsid interior will 
be simultaneously impeded. Release of the blockage by a further conformational change 
in the tail proteins then allows resumption of water flow and thus of DNA ejection. 
Interestingly, the apparent activation energy for releasing the block in the step-wise 
ejection process of T5 DNA ejection has been experimentally shown to be independent 
of the amount of DNA remaining in the capsid [371I70] . 

In summary, the hydrodynamic model is fully compatible with the experimental 
data for A and for SPPl DNA ejection in vitro. Like thermodynamics-based models 
(e.g., continuum mechanics or molecular dynamics), the hydrodynamic model requires 
an ad hoc assumption to explain the step-wise process of T5 DNA ejection. However, 
the hydrodynamic model is not dependent on the critical assumption — used in Pejection- 
based models — that the encapsidated DNA remains at equilibrium at all times during 
the ejection process. Furthermore, as we argue below, only the hydrodynamics model 
explains how complete genome ejection can be achieved in the face of an opposing force. 
This is the situation in natural infections of bacterial cells. 

4-2. DNA ejection in vivo 

Before discussing specific mechanisms of phage DNA ejection in vivo, it is instructive 
to consider some salient points about phage infections. 

4-2.1. Phage infection of bacteria Largely because of influential textbooks and, until 
recently, of only little experimental information, it is not widely appreciated that all 
phages eject proteins into the cell [69l|7l]. At a minimum, the protein(s) comprising the 
tail plug must be removed, and for long-tailed phages (like T4, T5, SPPl and A) the 
tapemeasure protein, which determines the precise length of the tail, must be ejected 
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in order to allow the phage genome to pass through the tail tube. Short tailed phages 
may eject internal proteins into the cell to extend — at least functionally — their tail 
so that it can span the infected cell envelope [57j . 

Some phages eject many different protein molecules from their capsid, some 
necessarily before, others perhaps after, genome ejection. To give just two examples, 
T7 virions ejects ~ 75 molecules representing five different protein species into the cell 
prior to DNA penetration of the cell cytoplasm [72], and T4 virions eject ~ 1000 IP 
(internal protein) molecules into the cell [73] . Clearly, any model purporting to explain 
even dsDNA phage genome ejection in vivo must also accommodate ejection of virion 
proteins. 

Phage virions can also eject more than one DNA molecule [7HI75] . However, as the 
diameter of the channel through the portal complex and tail tube (inner diameter ~ 40 
A) cannot accommodate more than a single DNA helix, the leading end of the second 
(and perhaps subsequent molecules) must find the exit channel without a built-in guiding 
vectorial force. Genome ejection by single-stranded DNA (ssDNA) or RNA viruses, 
whose nucleic acid is packaged at much lower densities than dsDNA phages, should also 
be explained by models describing phage DNA ejection. Packaged genomes of these 
phages are not thought to be under pressure, and their mode of genome ejection has 
therefore not been considered by continuum mechanics or molecular mechanics models. 

The E. coli cytoplasm has a positive osmotic pressure of several atm above the 
environment (under various growth conditions the pressure has been estimated to vary 
between 2 and 15 atm, with 3.5-5 atm being commonly accepted values [76l|77]). This 
osmotic pressure gradient is counter-balanced by a hydrostatic pressure differential 
(turgor) [68] that enables the cell to enlarge during growth ^S\- Gram-positive cells, 
such as Bacillus subtilis, the host for phage SPPl, have much higher turgor: ~ 19 
atm [79j. If turgor is assumed to provide an opposing force to phage genome ejection 
into the cytoplasm, then a Pgjection-based mechanism cannot possibly explain complete 
genome ejection into cells. This problem was recognized in the first osmotic suppression 
of genome ejection experiments in vitro [38], but has often been ignored in subsequent 
general statements about how phages infect cells in vivo. Interestingly, the osmotic 
pressure of the cytoplasm declines as bacteria enter stationary phase, which should allow 
better genome penetration by a Pgjection-based ejection mechanism. However, stationary 
phase bacteria are not readily infected by most phages [91], which infect exponentially 
growing cells, with their higher turgor pressure, much more efficiently. 

A bacterial cytoplasmic membrane is composed of a phospholipid bilayer that is 
impermeable to most molecules other than water and glycerol. An electrical potential 
is maintained across this membrane. A strong osmotic gradient also exists between 
the external medium and the cell cytoplasm. The higher internal pressure (turgor) is 
necessary for bacteria to enlarge and undergo cell division. In addition, at all times, a 
K"*" concentration gradient [K+], with [K+] high inside, is maintained in cells (a reverse 
gradient of Na"*" concentration also usually exists). These conditions are perturbed when 
a bacterium is infected by a phage as a direct connection between the cytoplasm and 
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the external medium is opened. This connection passes through the phage capsid and 
tail, and if it is open, and, as we have argued, water will flow from the external medium 
into the cell to neutralize the overall osmotic gradient. Simultaneously, K"*" will flow 
from the cytoplasm into the external medium. Furthermore, the membrane potential 
can no longer be maintained. 

Advocates of Pejection-based genome ejection in vivo have made various ad hoc 
suggestions to explain how an entire phage genome can enter the cell cytoplasm, 
including the involvement of cytoplasmic, non-sequence-speciflc, DNA-binding proteins 
or condensation of the phage DNA in the crowded cell cytoplasm. There is however, no 
direct experimental support for any of these ideas. Condensing the entering phage 
genome in the bacterial cytoplasm would effectively prevent its transcription, but 
efficient — and immediate — gene expression is precisely the primary strategy of 
phage infections. Furthermore, with the specific exception of second-step transfer of 
T5 DNA, it is hard to imagine how entering phage DNA could effectively compete 
with the 100-fold higher DNA concentration of the bacterial chromosome for those 
proteins in order to complete genome internalization in a kinetically reasonable time 
frame. T5 second-step transfer requires the prior synthesis of two T5-encoded DNA- 
binding proteins, and early T5-encoded nucleases completely degrade the bacterial 
chromosome [80], thereby removing competitor DNA. Nevertheless, sequence-specific 
DNA-binding proteins, which have no or fewer binding sites on bulk chromosomal DNA, 
and can therefore effectively bind to incoming phage DNA in vivo, have been shown to 
catalyze genome internalization of phage T7 and its relatives [69], and also of N4 [20] . 

4-2.2. The hydrodynamics model is consistent with the physiology of phage infection 
DNA ejection in vivo and in vitro are actually quite different. The major distinction 
is that when the tail is unplugged in vitro, the environment at the distal end of the 
tail tube is the same buffer that surrounds the capsid. In contrast, during a natural 
infection the distal end of the tail tube is in the cytoplasm of the infected cell, while the 
environment of the capsid is the growth medium. Consequently, the osmotic pressure 
of the medium that surrounds the phage capsid is different from that at the opening of 
the tail tip. 

We denote the osmotic pressure of the culture medium by vTout, and that of the 
cytoplasm of the infected cell by vTcyt; note that vTcyt > vTout as the bacteria need to 
maintain positive turgor in order to grow [77]. Although in Pejection-based model of 
phage genome ejection turgor is assumed to provide a resisting force against Pejection, 
it is the bacterial Achilles heel that actually promotes complete genome ejection by 
a hydrodynamics mechanism. When the tail plug is removed by interaction with its 
receptor, water will fiow up the osmotic gradient — one that exists between the growth 
medium and the cell cytoplasm — from the growth medium, through the capsid shell 
and tail tube into the cell cytoplasm. This directional water fiow can drag the DNA 
into the infected cell. This mechanism is schematically shown in Fig. [6l It should be 
noted that bacterial cells respond immediately to changes in Hcyt and they will always 
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maintain some turgor even during phage infection. 

When bacteria enter stationary phase, their turgor pressure drops considerably. 
This reduces the osmotic gradient between the external medium and the cytoplasm, in 
turn this reduces the strength of water flow from the culture medium into the cytoplasm. 
This may explain why bacteria in stationary phase are not easily infected by most 
phages. Notably, T7 is an exception, infecting stationary (and starved) cells at normal 
efficiency [91], even though subsequent phage development may be compromised by a 
lack of biosynthetic capacity in the cell. 



water entering through the capsid shell 




Figure 6. Schematic of DNA ejection in vivo, (a) Mature phage adsorption on a 
bacterial cell, (b) The plug blocking the tail is opened or removed as the cell envelope 
is breached. Water and small solutes enter the phage capsid from the culture medium 
and drains through the tail tube into the cytoplasm of the infected cell, dragging the 
DNA along with it. 

Following infection by all but one of the phages we are aware of that have been 
tested, there is a transient drop in membrane potential and a transient release of intra- 
cytoplasmic ions into the external medium [811182] . Both occur during the period that 
DNA is thought to be ejected from the virion into the cell cytoplasm. In order to 
evaluate the role of the membrane potential during T4 infection, an extensive series 
of electrochemical studies that monitored ion fluxes were conducted (Reviewed in 
Refs. [HnilHl]). Experiments were often combined with infective center assays after 
Hershey-Chase blending of phage-cell complexes, allowing estimates of 10^-10"^ bp/sec 
for the rate of T4 genome internalization in vivo [8T |[83] - [85] . Furthermore, T5 releases its 
DNA in two distinct steps pO]. During each step, cytoplasmic K"*" leaks into the external 
medium, while during the pause between the steps no leakage occurs [86]. Thus, ion 
leakage is associated with phage DNA ejection, and if ions are leaking out, water should 
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be flowing from tlie external medium through the phage virion into the cell cytoplasm 
to reduce the osmotic gradient. This is exactly what is expected if hydrodynamic forces 
are dragging the phage genome into the cell. 

There are conflicting published data on membrane depolarization and ion leakage 
for phage T7, but it is now clear that T7 infection does not result in either phenomenon 
[571 EH], indicating that no open channel that could allow water flow exists between 
the environment and the infected cell cytoplasm. These are significant observations, 
because, unlike other phages, the internalization of T7 and related phage genomes is 
entirely catalyzed by molecular motors. These motors are enzymes that utilize cellular 
energy to transport the entire genome into the cell, each functioning at a constant rate 
regardless of how much DNA remains in the phage head [56ll89] - [93] . Thus, in the absence 
of hydrodynamic forces, energy-requiring enzymes are necessary to effect phage genome 
internalization by the cell. 

The hydrodynamic model of DNA ejection does not require secondary ad hoc 
mechanisms to complete genome transfer into the cell. Because the water flow is 
determined by the values of vTcyt and vTout, the ejection process is (relatively) independent 
of DNA condensing agents. A bacterial cell responds immediately to changes in vTcyt and 
will maintain some turgor even during phage infection. Thus, hydrodynamic flow from 
the external environment through the phage particle into the cell cytoplasm will continue 
until genome internalization is complete and the channel is closed. (It is not known how 
this occurs in any phage system, but all models of phage DNA ejection must invoke this 
step in order to prevent a permanent loss of the cellular membrane potential. If the 
potential is totally collapsed, the cell would die and no phage progeny would ever result.) 
Furthermore, although low concentrations of DNA condensing agents have a major effect 
on TTin, they have little effect on either vTout or TTcyt (e.g., intracellular levels of polyamines, 
which are critical for ribosome and chromosome stability are tightly regulated). This 
means that phage DNA ejection in vivo should be relatively independent of the presence 
of DNA condensing agents in the growth medium. Indeed, ongoing experiments in the 
laboratory of one author (IJM) suggest that the latent periods and burst sizes of various 
phages (including parallel experiments with A and the deletion mutant Ab221, the latter 
having a substantially reduced value, relative to wild-type A, of vTin ) are not affected by 
the presence of 1 mM spermine in the bacterial growth medium. 

The hydrodynamic model also provides a mechanism for the ejection of single- 
stranded genomes and proteins into the cell during infection. This is important as all 
phages eject protein molecules into the infected cell, and a significant fraction of phages 
do not contain dsDNA packaged at the high density of ~ 500 mg/ml. Any molecule in 
the path of water flowing along the osmotic gradient between the external medium and 
the cell cytoplasm can be driven down the tail tube and into the cell. As this water 
flow is necessarily vectorial, nucleic acids or proteins to be ejected do not need to be 
positioned in the tail tube in the mature phage particle. Neither do they need to possess 
internal energy to drive their ejection from the capsid. 

In summary, with the caveat that, to explain the multi-step process of T5 
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and perhaps other phages (029 is one known example [SS]), an additional ad hoc 
mechanism of protein conformational changes that temporarily block water flow from 
the external fluid into the cell cytoplasm, in turn temporarily stopping DNA ejection, 
the hydrodynamic model of phage DNA ejection is consistent with — and importantly 
can explain — all the observations we are aware of that have been made using any 
phage-host combination in vivo. 

However, theories need to be tested experimentally. Determining the kinetics of 
phage genome internalization into infected cells would go a long way to supporting or 
refuting competing models. The internal capsid pressure models predict that the late 
stages of genome ejection should slow towards a zero rate as the residual pressure equals 
that in the cell cytoplasm, and that a secondary process using, for example, polyamines 
or DNA-binding proteins to interact with the entering the DNA completes the process. 
A secondary process is likely to occur with different kinetic parameters and may be 
experimentally detectable. Conversely, the hydrodynamic model would predict a largely 
constant rate of genome internalization, as the osmotic pressure gradient between the 
cell cytoplasm and the environment is unlikely to change substantially while viral DNA 
enters the cell. Unfortunately, it has proved very hard to measure directly the in vivo 
kinetics of genome entry outside the T7 family of phages, although some measurements 
have been made with phage N4 [96], and experiments with both A and T5 are being 
initiated in the Molineux laboratory. However, it may be possible, using a combination 
of polyamines and external osmolytes, to increase the cytoplasmic osmotic pressure 
above that inside the mature phage head. If phage infection still occurs normally under 
those conditions, it would then provide irrefutable evidence that internal capsid pressure 
is not important for infection, although it would not provide any direct evidence for the 
hydrodynamic model. 
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